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Abstract. One of the challenges of the Swift era has been accurately determining Epf,ak for the 
prompt GRB emission. RHESSI, which is sensitive from 30 keV to 17 MeV, can extend spectral 
coverage above the ^wZ/f-BAT bandpass. Using the public Swift data, we present results of joint 
spectral fits for 26 bursts co-observed by RHESSI and Swift-QP^ through May 2007. We compare 
these fits to estimates of Epgak which rely on BAT data alone. A Bayesian Epgak estimator gives 
better correspondence with our measured results than an estimator relying on correlations with the 
Swift power law indices. 
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GRB PROMPT EMISSION SPECTROSCOPY WITH RHESSI 

The Reuven Ramaty High-Energy Solar Spectroscopic Imager (RHESSI) fl'] is a dedi- 
cated solar observatory. RHESSI's nine germanium detectors are sensitive from 30 keV 
to 17 MeV, with excellent resolution in energy (1-5 keV) and time (1 binary jUs) 0]. Each 
of the nine coaxial detectors is electronically segmented into front and rear segments. 
Because the detectors are unshielded, RHESSI observes emission from astrophysical 
sources like GRBs with a ~2;r field of view. 

We perform Monte Carlo simulations using the MGEANT package Jsl] to determine 
RHESSI's response to off-axis sources like GRBs. RHESSI's response varies with off 
axis angle, so we create responses every 15 degrees. For each response, we simulate 
monoenergetic photons in 192 logarithmic energy bins ranging from 30 keV - 30 MeV. 
Since RHESSI's per-detector response also varies during the spacecraft's four second 
spin period, we bin the annular response in six azimuthal bins and weight these bins by 
the total burst lightcurve to create the final response. 

The RHESSI data are extracted in SSW-IDL. We fit and subtract a time-varying 
background, allowing forpossible periodic modulation with the spin period. We perform 
spectral fitting in ISIS [4], a forward-fitting package analogous to XSPECQ, which is 
extensively programmable and allows computation of rigorous fluence error estimates 
via exploration of the parameter space. 

Energetic charged particles over time have caused radiation damage to RHESSI's 



http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/ 



germanium detectors. This radiation damage causes broadened spectral lines due to hole 
trapping and a general loss of active volume. In this work, we restrict our analysis to 
those detectors which do not exhibit signs of radiation damage. At the time of Swiff s 
launch, six RHESSI segments were usable for spectroscopy; by May 2007 only two 
segments remained undamaged. In November 2007, RHESSI underwent an annealing 
procedure to reverse the effects of the radiation damage. The anneal restored some lost 
sensitivity, but analysis of future bursts will require more sophisticated modeling of the 
remaining effects of radation damage on RHESSI's spectral response. 



RHESSI-BAT JOINT FITS 

We attempted simultaneous spectral fitting for all RHESSTobserved GRBs appearing in 
the first BAT Catalog [5]. Of 46 candidate bursts, 26 had sufficient RHESSI counts for 
spectral analysis and produced acceptable joint fits. We selected analysis time intervals 
manually from the RHESSI lightcurve using a S/N criterion. The resulting intervals 
were usually shorter than those used in the BAT Catalog. We generated BAT spectra and 
responses for our intervals with the standard analysis procedures H 

Typically, the joint fits did not require a normalization offset between the RHESSI 
and BAT data. (Of the four bursts that needed an offset for acceptable fitting, one was in 
a period of highly modulated RHESSI background, and three were after December 2006 
when radiation damage was becoming severe.) For the RHESSI data, we generally fit 
over the full 30 keV-17 MeV energy band. For bursts coming from the rear of RHESSI, 
we raised the lower energy bound to ~60 keV, as the additional passive material of the 
RHESSI cryostat can influence the low energy data. For GRB 061007, we omitted the 
RHESSI data above 3 MeV; there were no significant counts above that level, but the 
residuals showed systematic deviation which biased the fit. 

Our results show good correspondence with comparable fits reported by Konus-Wind 
and Suzaku-WAM, which also are sensitive in the MeV range [e.g.,laLZil8i|9i[lO|]. For 



16 of the 26 bursts, the joint RHESSI-BAT best fit found additional model parameters 
(Epeak and/or /3) compared to the BAT-only fit. 
We report the results of the joint fits in Table[T] 



TESTING Epeak ESTIMATORS 

The peak energy Epeak of the prompt GRB spectrum is crucial to determining overall 
burst energetics, and it plays a key role in several proposed luminosity indicators. 
The narrow passband of Swift-BPJ' prevents determination of Epg^k for many bursts. 
Accordingly, a number of attempts have been made to infer Ep^ak from the BAT data 
alone. Butler et al. [11] used a Bayesian fit method with priors determined from the 
BATSE catalog to estimate Epeak- Zhang et al. [12] derived an Epeak - BAT power law 
index correlation using hardness ratios (see also [|l3n ). In Figure [H we compare the 



http://swift.gsfc.nasa.gov/docs/swift/analysis/threads/bat_threads.html 
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FIGURE 1. Compai-ison of the measured RHESSI+BAT 

Epeak values with those predicted by ll ill 
(left) and those found with the correlation of IJ2.1 (right). Points marked with a cross have fit j3 > —2, and 
hence Epg^^i^ is only a formal value. Points in the right plot marked with a diamond have BAT-only power 
law indices outside the range (—2.3 < a < — 1.2) likely to yield accurate predictions of Ep^ak 1131]. The 
overplot lines represent equality of the measured and predicted values. 



predictions of these models to our joint fit results. 

The Bayes model shows good correspondence with the measured values. Above ~600 
keV, the predicted values of Epeak tend to be low, although their error bars reach near the 
measured values. This deviation is exaggerated somewhat, as our data are for shorter, 
more intense (and typically harder) burst intervals than used in [1 1]. 

Comparison of the Zhang et al. correlation-predicted Epgak to our measured RHESST 
Swift values shows that this correlation appears to systematically underpredict the mea- 
sured value of Epgak, especially at high energy. 
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TABLE 1. Best joint fit parameters. Spectral models considered are a Band function Ill4ll . cutoff power law 
(CPL), and simple power law (PL). Errors are quoted at the 90% confidence level. Fluence is in the 15 keV-10 
MeV band, and analysis times are relative to the BAT trigger time. 
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